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Here we try and delienate the properties of the function that corresponds to fluctuations 
in the momentum distribution. The quantity denoted by N(k, k ) is quite an interesting 
", object. It satisfies various elegant sum rules and is also quite singular in some respects. All 

these properties are brought out and a formal connection is found between this object and 
the momentum distribution of the interacting and non-interacting many-fermion systems. 
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' This exercise is quite general in that it does not refer to any particular hamiltonian. It is 
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also quite useful since in an earlier preprint(cond-mat/9810043) we showed how to compute 



the spectral function and single-particle lifetime of homogeneous Fermi systems where the 
only undetermined quantity was this function N(k, k ). 

I. NUMBER FLUCTUATIONS 

Without further ado, let us define the object of interest, and try and enumerate its many interesting 
properties. The number- fluctuation N(k, k ) is defined as, 

iV(k,k') = Kn k O-(n k )KO (1) 
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here n k — c k Ck- Ck and c k are Fermi fields. It is useful at the outset to point out some rather trivial facts 
q | about this object, 



X>(k,k') = (2) 

k 

This has to do with the fact that the total number of particles is conserved. But for those who are skeptical a 
simple demonstration may settle any lingering doubts. Also this approach will be used for illustrating more 
substantive issues later on, hence its introduction at this stage is quite appropriate. 

J2N(^') = (Nn u ,)-(N)M (3) 

k 

Let {|*}o} be the eigenstates of the noninteracting Fermi system with total number of particles equal to No. 
And {\i)} be the eigenstates of the interacting Fermi system with total number of particles equal to Nq. We 
may write, 



<^M = |E e ~^i^vb> (4) 
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here Z = J^j e~^ Ej is the canonical partition function. But, 
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(J\N rv \j) = Y,V\N \i)(i\n k , \j) = N (j\n^ \j) (5) 

i 

Since, (N) — N also, the result follows. A similar result holds for the "first moment" namely, ^ k k/V(k, k ). 
Since P = ^ k k«k is the total momentum of electrons, the same argument gives us, ^kk^(k, k ) = 
0. But what about those cases when the eigenstates of the interacting system are not eigenstates of the 
operator P ? This case for example corresponds to the situation when translational invariance is broken 
spontaneously (Wigner crystal). The answer is as follows. 

(j|P rv|j) = E01 p Ko^olrvli) = £n»,P i0 |<io|j>| 2 (6) 

io io 

here {\i0)} are the states of the free theory. Therefore Pjo = for all zO. But then the overlap (io I J ) may 
be singular. By that I mean it may be actually infinite over a set of measure zero. Thereby orthonormality 
is not compromised and yet J] k kZV(k, k ) may fail to vanish. Let us prove another interesting property of 
this function. 

7V(k,k) = K)-(n k ) 2 = (n k )(l-(n k )) (7) 

A quick proof of this involves making the following observation, n k = ra k . But this equation has to be 
interpreted with care. It is not an operator identity in the usual sense. If it were, many strange (and 
incorrect) results would ensue. For examplc( |q| << |k|), 

(«k+ q ) « (rik> + ((q-Vk)rik) 

= (n k ) + <(q.V k )n£) = (n k ) + 2<n k (q.V k )n k ) 
(because (q.V k )n k w n k+q - n k and [n k , (q.V k )n k ] = 0) 

= (n k ) + 2(n k (n k+q - n k )) = 2(n k n k+q ) - (n k ) 

This means, 

(^k^k+q) - (nk)(n k + q ) = -^(-(n k+q ) - (ra k ) + 2(n k )(n k+q )) 

We can see quite easily that this equation is false by applying it to the noninteracting system at zero 
temperature. The left side is identically zero and right side may not be. However, we may still apply the 
ideas introduced earlier to write down some useful formulas for these objects(N's). Then in what sense is 
n k = n k true ? The answer is, in the sense of matrix elements. 

(io\K\jo) = X^°l nk l fc °)( fc o|rik|:/o} = n k °n k ^5i ,fc 5j ,fc = {K) 2 6 io , jo (8) 

ko ko 

Since the c-numbcr is either zero or one, the result follows. A further subtle issue comes to mind. In the 
case of noninteracting systems at finite temperature, iV(k, k) ^ but N(k, k ) = for k ^ k . Therefore 
N(k, k ) is discontinuous at k = k . In particular, 
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L?m |qUq|'-o^( k - q, k - q ) ^ JV(k, k) (9) 
Having done all this we now proceed to write down formulas for the off-diagonal elements of N. 

N(k,k) = i£e-^<i|n k rv|i) - (n k )<n k (10) 

i 

here \i) are the eigenstatcs of the full theory. We may now decompose 

(i\n k n k ,\i) = "£{i\n k \j ){j \n k/ \i) =Y,ni°n$\(j \i)\ 2 (11) 



Further, 



^e-p*®^) = ^E^E^KWOoi^') = ^E e ~^ 4 E4°iwi 2 ( 12 ) 
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Let us introduce the inverse of the "matrix" {n^ } 

J2<°( n '% = 5 ^' ( 13 ) 



jo 



Therefore, 



E<°( n " 1 )iO=W (14) 



J2( n *)( n ~% = -E^i^i 2 ( 15 ) 



(n k TV> = D|E e_/JBl IO'o|i>| a K°ni' ( 16 ) 



Therefore, 



<»k«k'> ^E^H^'IKX" ( 17 ) 

Jo P 



Thus, 



iV(k,k') = ^EWt"" 1 )^ ~ K)K'> (18) 

jo P 

The above formula is exclusively in terms of the (nonideal) average (n p ). The rest of the quantities are known 
except the "inverse" of n 3 ,, may be hard to compute. This preprint together with our earlier one (cond- 



mat/9810043| ) presents a viable alternative (hopefully) to other approaches at computing single-particle 



properties of Fermi systems. 
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